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Abstract. We investigate the contribution to the angular correlation coefficients of 
the neutron beta decay within the R-parity violating (RPV) minimal supersymmetric 
standard model (MSSM). The RPV effects contribute to the scalar interaction at 
the tree level. The coefficient of the effective scalar interaction of the nucleon is 
determined in terms of the nucleon mass difference. On the basis of the recent update 
of the analyses of the superallowed Fermi transitions and the recent measurement of 
transverse polarization of the emitted electrons at PSI, we deduce new upper limits 
on the RPV couplings and compare them with those obtained through pion decay 
and atomic electric dipole moment (EDM). We also give a comprehensive analysis of 
angular correlations sensitive to the RPV interactions. 
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1. Introduction 

In recent years, the search for new physics "Beyond the Standard Model" has been 
performed in many ways. A particular attention is paid to the low energy approach, 
which consists of observing small discrepancies between low energy phenomena and the 
Standard Model (SM) [H Ej [3l HJ |5] . Among several others, neutron decay experiments 
are expected to be very promising for new physics search and are actively planned or 
under operation by many groups at PSI, LANSCE, ILL, J-PARC and so forth. 

On the theoretical side, the Minimal Supersymmetric Standard Model (MSSM) 
[6] is known to be the leading candidate of the new physics beyond the SM. We must 
impose on the MSSM a new parity called R-parity (Rp = (_ i)3B+i+2s^ ^.^ conserve the 
baryon and lepton numbers, in contrast to the SM in which these quantum numbers 
are automatically conserved due to the gauge invariance. The violation of the R-parity 
causes many exotic phenomena such as the proton decay. Although its conservation 
is generally assumed, this assumption is completely ad hoc. The RPV interactions 
may exist which could be unveiled by upcoming high-precision experiments. Until now, 
many allowed regions and upper bounds have been established for single and combined 
parameters of the RPV interactions by means of the weak processes [TJ [HI El fTO] . 

Herczeg [1] analyzed the d — )■ ue'Ve transition by including the RPV interactions 
and pointed out that they generate two types of new interactions. One is the usual 
V-A interactions and the other is of the scalar and pseudoscalar type. It was shown 
that rather severe upper bounds on the RPV interactions [H] were obtained from the 
contribution of the pseudoscalar interaction of the pion decay r(7r — ei^e)/r(7r — )■ yuz/^) 
[T2] . The bound on the scalar interaction was also derived by the use of electroweak 
radiative corrections [I3] and the atomic EDM data [TTj . 

Recently, new data sensitive to the scalar interaction have been reported. Hardy 
and Towner [H] updated their analysis on the superallowed Fermi transitions and gave 
an improved limit on the real part of the scalar interaction between the hadron and 
lepton currents. Meanwhile, the experiment at PSI [15] has measured the transverse 
polarization of electrons emitted in the neutron beta decay, and has given observables 
sensitive to both the real and imaginary parts of the scalar interaction. 

As the scalar interaction is sensitive to the RPVMSSM contributions at the tree 
level [1], the present analysis of neutron beta decay will provide useful constraints on 
RPV interactions complementary to those obtained from the pseudoscalar interaction 
of the pion decay and the atomic EDM via one-loop corrections. 

The purpose of the present paper is to perform the quantitative analysis of the RPV 
contributions to the neutron beta decay in the light of [11] and [I5]. In order to explore 
a possibility to investigate the scalar interaction from the neutron beta decay, we use 
the formula of angular correlation coefficients given in literatures [T6l [T7] . We notice 
here assumptions we made in this paper, which are the foUowings: (a) The R-parity 
conserving sfermion sector of the MSSM is assumed to be flavour diagonal and with no 
CP phases, (b) The RPV sector of the MSSM does not contain any soft SUSY breaking 
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terms. 

This paper is structured as follows. In sec. |2] we first construct the effective 
interaction of the neutron beta decay within RPVMSSM. In sec. [3], we explore 
the angular correlation coefficients by including neutrino momentum in the electron 
polarization term. In sec. IH we obtain new constraints on the RPV couplings from the 
recent data of nuclear beta decay. Sec. [5] is devoted to summary. 



2. Effective interaction from R-parity violation 

2.1. R-parity violating Lagrangian 

The first step of the estimation is to construct the tree level amplitude of the quark 
beta decay. The relevant interactions can be obtained from the RPV part of the 
superpotential which can be written as follows: 

with i,i,k = 1,2,3 indicating the generation, a,b = 1,2 the SU{2) indices, l,m,n = 
1,2,3 the colour indices. L and E'^ denote the lepton doublet and singlet left-chiral 
superfields. Q, and D'' denote the quark doublet, up quark singlet and down quark 
singlet left-chiral superfields, respectively. The third term in ([T]) is baryon number 
violating, so is not relevant to our subsequent discussions. We therefore disregard it 
hereafter. By taking the F-terms of W-^, we obtain lepton number violating Yukawa 
interactions and also scalar four-point interactions. In our discussion, the scalar four- 
point interactions are irrelevant since we cannot construct diagrams which involve them 
at the tree level. The RPV interactions contributing to the neutron beta decay is as 
follows: 



ijk 

- e^^^.e'iPLVj - CLiCkPLVj - VjCkPhei^ 
~ Yl Kjk{d'^Rk^iPLdj + dLjdkPh^i + VidkPhdj 

ijk 

- i^i.e'lPLUj - eLidkPiUj - ULjdkPiei] 



h.c. 



(2) 



withP,. = i(l-75). 



2.2. Feynman diagrams 

We can now build the amplitude of the quark beta decay within RPVMSSM at the tree 
level. The Feynman diagrams are shown in Fig. [TJ 

Fig. [1] (a) is the SM contribution at the lowest order. Its expression is as follows: 

G 

MsM = -^Vud U-f^'il - -f5)d ■ 67^(1 - 75)z^e , (3) 
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(a) (b) (c) 

Figure 1. Tree level amplitude of the quark beta decay within RPVMSSM. 



where Gp = (1.16637 ± 0.00001) x 10-^GeV"^ [H] is the Fermi constant and Kd = 
0.97418 ± 0.00027 [19] is the CKM matrix element. 

In the tree approximation, the RPV contributions to the d ^ u + e~ + process 
are the selectron(ei) (Fig. 1 (b)) and the squark((iR) (Fig. 1 (c)) exchange mechanisms 
[1]. The amphtude of the selectron exchange between the quark and lepton in Fig. [1] 
(b) can be written as 



(Alii — Ajii)A^^^ _ , , \ 7 -/I \ 

Me^ = 2^ —2 —U (1 + 75) ■ e (1 - 75) I/e 

i=2,3 ^"^eL^ 

= E + 75)rf ■ e(l - 75)^^e , (4) 

i=2,3 ^^eL, 

where we have used the antisymmetry of the couphngs Xijk under the interchange of the 
first and second index {Xijk = —Xjik)- It is interesting to note that A^ez, is of the type 
of the scalar interaction. 

The RPV couplings are constrained by extensive analyses on current experimental 
data. The upper limits of the RPV couplings used in (jlj) are given by the first four entries 
in Table [TJ The upper limits of Aiii, X'^ for z = 2, 3 are obtained from charged current 
(CC) universality and decay ratios r(r — ?■ ez/P)/r(r — )■ fiiyi^), r(vr — ei/)/r(7r fiu) 
and r(r — 7rz/^)/r(7r — )> fiu) [TJ [H [9]. As we see, the magnitude of A^e^ is small 
compared with that of A^sm. This scalar interaction is, however, potentially important 
since it could give rise to observables in neutron beta decay to which V — A interactions 
do not contribute. 

The squark exchange amplitude in Fig. [T] (c) can be written as 



-^..= E l^^:(l-75M-n(l + 75)e'^ 



=1,2,3 



E ^^7ni-75)rf-^e^7M(l + 75)e^ 



i=l,2,3 ■^^dn. 



E ^-^Ur{l-l5)d-n,{l-l>e. (5) 



i=l,2,3 ■^^dR.^ 



To obtain the final result in (E]), we have used the Fierz transformation. has the 
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Table 1. Current upper limits of (the magnitude of) RPV couplings [H [20] , 
Expressions with square brackets correspond to their masses in unit of 100 GeV. q 
and g refer to the squark and gluino, respectively. 



Coupling constants Upper bounds 




< 0.049 [mg^] 

< 0.062 [mg^i 




< 1.3 X 10-4. [m5]2[mg]i/2 

< 0.021 [m-s^] 

< 0.021 K^] 



V —A type interaction like the SM. Here the constraints on the RPV couphngs appearing 
in (E]) are given by the last three entries in Tabled! The couplings A^^j^ [i = 1,2,3) are 
constrained from the analysis of the double beta decay [20] and CC universality [3 13 [9] . 
In our analysis, we use the empirical value of Vud and axial vector coupling constant qa, 
therefore the effect of ([5]) is absorbed into those coupling constants. 

2.3. Effective interaction 

The next step is to build the effective interaction of the neutron beta decay process. To 
obtain it, we must construct the nucleon matrix elements from the quark beta decay 
amplitudes. The effective interaction can be written as 



where the first term is the SM contribution and the second term the contribution from 
Fig. [1] (b). The nucleon matrix elements are described by form factors as follows: 



where we have explicitly written the nucleon fields as n and p, and q is the momentum 
transfer. The form factor of the pseudoscalar interaction need not be introduced since 
its contribution vanishes in the non-relativistic limit. The induced form factors were 
also neglected. In this paper, we use the constants gy = fl'y(O), = fl'A(O), gs = dsi^) 
since q"^ is much smaller than the nucleon mass. Here recall that gv = 1 from the CVC 
and gA = 1-2739 ± 0.0019 from [2l|. 

The determination of gs needs more involved discussion. Within approximate 
isospin symmetry, the matrix element {p\ud\n) = {p\uu—dd\p) . According to Gasser and 
Leutwyler [22], the isospin asymmetric part of the nucleon form factor can be written 
as 



Hp = HsM + 



(6) 



{p\u'y^d\n) = gv{q^)pi^n , 
{p\u-i^l->d\n) = gA{,q^)Pl^lbn , 
{p\ud\n) = gs{q^)pn , 



(7) 
(8) 
(9) 



5 



rUfi — m. 



'U 



((p\uu — dd\p) , 



(10) 



2Mjv 
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where M^r indicates the nucleoli mass. To the lowest order expansion in up and down 
quark mass difference, the neutron and proton mass difference equals 5: 

S = Mn-Mp = (2.05 ± 0.30) MeV , (11) 

where the mass difference Mn — Mp are of pure QCD origin, i.e., without electromagnetic 
contribution. From (191), (fTOl) and (fTTl) we obtain 



gs = — ^ = 0.49 ±0.17, 12 

ma - TTLu 

where we have used = 9.3 ± 0.9 MeV, mu/md = 0.553 ± 0.043 [23]. The error in 
f ll2p is mostly due to the ambiguities of the current quark mass. Note that the values of 
Qs obtained previously by Adler et a/.[2l], and by Wakamatsu [25] are consistent with 

m- 

We can now construct the effective interaction of the neutron beta decay within 
RPVMSSM at the tree level, which is given by 

Gf - 

HsM = —VudP1^{gv - 9Alb)n ■ 67^(1 - 75)z/e (13) 

v2 

and 

i/^=Cspn-e(l-75)i/e + h.c. (14) 

where 

3. Angular correlations of the neutron beta decay 

The formula for the angular correlation coefficients of the beta decay by taking into 
account of the polarized neutron and electron is given by Jackson, Treiman and Wyld 
(JTW) [16], and extended by Ebel and Feldman [17] as follows: 

u{Ee, l^e, f^,.) ocl + a^p^ + 

. ^ j ,Pe . r^Pu . ^PeX Pu] 
+ (Jn ■ < A— + B— + U- 



Ee Ey E(,Ey J 
I Ee Ey Ef. + nif, Ef,Ey E^E^ 



[ Ee + m Ee Ee ] 

[ J^e-t^u i^ei^u i^e^p 

Ee + me E,E, J ^ ' 

Here Pe and p^ are respectively the momenta of the electron and the neutrino, a„ 
the neutron spin polarization and the final state electron polarization. The angular 
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correlations of the last line (terms with S", T, U , V and W coefficients) are correlations 
with momentum of emitted anti-neutrino, polarizations of initial neutron and emitted 
electron, given by Ebel and Feldman [T7] . 

The angular correlation coefficients are evaluated using the main term of vector and 
axial vector current of the SM in ( fT3l) and the scalar interaction in (H^ll . Here we take 
static approximation of nucleon current and keep only the leading order contributions. 
The results are summarized in Table [2] where we used 

A = ^, (17) 
9v 



= Gp,r , . o^o^ ^ (18) 



^Kd^7y(l + 3A2 
2Im{Cs) 



ai = ^ -. (19) 

^Kd^7y(l + 3A2) ^ ' 

There are quite a few terms which are sensitive to the scalar interaction (see Table 
12]). The coefficients fe, 5, U and i?, L, V, W probe respectively the real and the imaginary 
part of the scalar interaction. Among them, S", T, V and W , in addition to the 6, R and 
L which already studied by JTW, are particularly interesting observables. Since the 
main term of the SM (or V — A) does not contribute to those observables, there is a 
possibility to probe new physics through the scalar interaction. The effects of the final 
state interaction (FSI) have been in [161 l26l [27] . 

It is also noticed that the B correlation can be used to investigate the scalar 
interaction. Though V — A interactions contribute to one may be able to extract 
a small effect of the scalar interaction by using the extra energy dependence ^ in the 
scalar contribution [1]. 

4. Analyses 

As we have seen, the RPV interaction of ([2]) generates (|1]) which turns out to be scalar 
interactions of the neutron beta decay (fT4)) . Combining the current upper limits of Table 
[H we have the following bounds for each combination of RPV couplings in (fT5]) : 

|Ai2iA';nl < 0.0029K-JK-J , (20) 

|Ai3iA'3\i| < 0.0068K-JK-J , (21) 

where [msusv] = ''7^susY/(100 GeV). f l2U]) and fl^ will be compared with the analysis 
confronted with [13] and |15] . 

From the recent update of the analysis on 20 superallowed Fermi transitions. Hardy 
and Towner [T3] have given a new bound on the Fierz interference term bp which 
corresponds to b term in ( IT6]) . This one is a probe of the real part of the scalar 
interaction. By fitting Ft values of the superallowed Fermi transitions, they obtained 

bp Re{Cs) 



2 a 



-0.0011 ± 0.0013 , (22) 



V 
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Table 2. Angular correlation coefficients of the neutron beta decay 



Coefficients 


SM 


RPV 


a 


[I- A + 6A ) 


n 
U 





U 




A 

/± 


ZA\L — A) / \i- -T oA ) 


U 


TD 
ID 


ZA[L + A) / [1 + 6A ) 


Aan\me/ -tje) 


U 


U 


U 


n 


1 

— 1 



u 


TT 
11 




— CUR 


1\ 


(A - + ^A j 


CUR 


L 





ai 


N 


-(me/i;e)2A(l-A)/(l + 3A2) 


-XaR 


Q 


-2A(1-A)/(1 + 3A2) 


XaR 


R 





— Xai 


S 





XaR 


T 


-2A(1 + A)/(1 + 3A2) 





U 





-XaR 


V 





—Xai 


w 





Xaj 



where Cy = Vud^dv- From this, we obtain the following limit on the combined RPV 
couplings: 

Re( Y: AiaAfn) = (7.2 ± 8.5) x lO"^!^ , (23) 
\=i^3 ^ [9s] 

where [gs] = 

The complementary information on the imaginary part of the RPV couplings can 
be obtained from the R correlation of the neutron beta decay distribution. The angular 
correlation R, which is the triple product of the polarization and the momentum of 
emitted electron, and polarization of the initial neutron, has been recently measured at 
PSI [15] as 

i? = 0.008 ±0.011 ±0.005. (24) 
The contribution of FSI to R can be estimated as 

Rfsi = 8.57 X 10"^ X — , (25) 

Pe 

according to the formula of JTW [16]. Here rrtg and Pe are respectively the mass and 
the magnitude of the spatial momentum of the emitted electron. The effect of FSI is by 
an order of magnitude smaller than the experimental data f l2^ . Using the experimental 
value f l2^ neglecting FSI, we obtained the following constraint on the imaginary part 
of the coupling constants, 

Im I XiiiK*u] = (-0.012 ± 0.017 ±0.008)^^^^^ . (26) 



vi=2.3 



[9s] 
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Im(Ei=2,3 ^lil-^'iil) 



-0.025 



-1 X 10 




2 X 10 



0.05 



Figure 2. Allowed region of ^ 
gs = 0.49. 




We have set mg, = 100 GcV and 



The obtained limits on the RPV couphngs in this analysis are summarized in Fig. 2. 
Here we have taken all s-particle masses to 100 GeV. As we can see, a strong constraint 
of the real part of the RPV parameters is obtained from b. On the other hand, the 
constraint of the imaginary part of the RPV parameters might be comparable to (|20ll 
and (12T]) . However, the limits obtained in ( l23l) and ( l26l) depend on the mass mg^ while 
the current ones ( f l20|) and fl2T|) ) depends on the masses mg^j and rn^^. Therefore (l23l) 
and (l26i) give a new type of constraint. With the progress of experiment on i?, one will 
be able to further constrain the imaginary part of the coupling constant. 

It is of interest to compare our constraints, ( !23l) and ( |26l) with those obtained 
previously. It has been pointed out in [13] that usual electroweak radiative corrections 
could transform the T-odd, P-even interactions into a T-odd, P-odd ones, and that one 
could put limits on the T-odd, P-even interactions by using the rich data on the part of 
T-odd, P-odd interactions. Herczeg [11] considered such radiative corrections in RPV 
models and confronted the induced interactions with EDM data. The analysis depends 
on the sparticle masses, but in general more severe constraints can be made available 
than ([26]). 

It has been pointed out in [12] that the ratio r(7r — )■ ev)/T{Ti — t- /iz/) is a sensitive 
proble to the pseudoscalar type du — ?■ ev^ interactions. For the case of the RPV 
contribution The scalar and the pseudoscalar interactions are of the same strength, 
and one can deduce an upper bound on Cs/gs by employing the data of the ratio. The 
upper bound obtained in such a manner is a fraction of the value ( 123|) . 
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5. Summary 

To summarize, we have studied the contribution of the RPVMSSM to the angular 
correlation coefficients and the Fierz interference term in the neutron beta decay. We 
have deduced new bounds on the real and imaginary parts of the combinations of RPV 
couplings as in ( l23l) and ( l26l) through the scalar interaction ( 1141) of the beta decay. These 
are based on the recent data of [H] and [I5] respectively. The obtained upper bound 
for Re{\iji\'*]^i) is consistent with a stronger constraint obtained from the analysis of 
pion decay [Him [12]. With the hypothesis of the single coupling dominance, f l23l) and 
( l26l) can be replaced by a constraint for either Ai2iA2\i or Ai3iA'3\]^. 

We have also discussed angular correlations S*, T, f/, V and W ^ which involve 
emitted electron polarization and neutrino momentum. In addition to the known ones, 
coefficients 6, 5, U and i?, L, V, W are respectively sensitive to the real and imaginary 
parts of the scalar interaction. Their precise measurements are therefore of great interest 
for further limiting the RPV couplings. Regarding the CP violating (imaginary) ones, 
the precise evaluation of the FSI contribution is also needed, as done in [26] or [27] . 
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